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Abstract. We describe a simple framework to assess the LISA scientific 
performance (more specifically, its sensitivity and expected parameter-estimation 
precision for prescribed gravitational- wave signals) under the assumption of failure 
of one or two inter-spacecraft laser measurements (links) and of one to four intra- 
spacecraft laser measurements. We apply the framework to the simple case of 
measuring the LISA sensitivity to monochromatic circular binaries, and the LISA 
parameter-estimation precision for the gravitational-wave polarization angle of 
these systems. Compared to the six-link baseline configuration, the five-link 
case is characterized by a small loss in signal-to-noise ratio (SNR) in the high- 
frequency section of the LISA band; the four-link case shows a reduction by a 
factor of \/2 at low frequencies, and by up to ^ 2 at high frequencies. The 
uncertainty in the estimate of polarization, as computed in the Fisher-matrix 
formalism, also worsens (but only for sources above 1 mHz) when moving from 
six to five, and then to four links: this can be explained entirely by the reduced 
SNR available in those configurations. In addition, we prove (for generic signals) 
that the SNR and Fisher matrix are invariant with respect to the choice of a 
basis of TDI observables; rather, they depend only on which inter-spacecraft and 
intra-spacecraft measurements are available. 

PACS numbers: 04.80.Nn, 95.55.Ym 

1. Introduction 

LISA (the Laser Interferometer Space Antenna) is a deep-space mission planned jointly 
by the National Aeronautics and Space Administration and the European Space 
Agency. LISA seeks to detect and study gravitational waves (GWs) in the mHz 
frequency band. It consists of three spacecraft flying in a triangular formation, whose 
relative positions will be monitored by way of laser interferometry [T]. In contrast 
to ground-based interferometric GW detectors, LISA will have multiple readouts, for 
the six laser links between the spacecraft (one in each direction across each arm). 
These data streams, properly time shifted and linearly combined, provide observables 
that are insensitive to laser-frequency fluctuations and optical-bench motions, and 
that have different couplings to GWs and to the remaining system noises [5]. This 
technique is known as Time-Delay Interferometry (TDI) . Early on it was realized that 
different TDI observables can be built from different subsets of the six inter-spacecraft 
laser measurements |3] , thus providing failure resistance against the loss of one or more 
measurements; furthermore, two or three linearly independent TDI observables can 
be used together to increase SNR and improve GW-parameter determination (this 
insight goes back at least to ^). 
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In this paper we pull together and complete the results scattered throughout the 
TDI literature (see references in [5]) to build a simple framework that can assess the 
multi-observable LISA science performance when all measurements are available, and 
in failure scenarios where up to two of the six LISA inter-spacecraft measurements 
and up to four of the six intra-spacecraft measurements are lost. (We do not analyse 
the case of three lost inter-spacecraft laser-links: in this scenario, LISA cannot make 
wide-band measurements of gravitational radiation, but narrow-band measurements 
are still possible at frequencies equal to integer multiples of the inverse round-trip 
light time between the two spacecraft that preserve bidirectional laser measurements 
0.) 

In section [2] we introduce our notation and our model of the LISA measurement; 
in [3l we describe a simple linear-algebra procedure to obtain the TDI observables 
that can be constructed in each LISA configuration; in Hj we prove that any choice 
of a basis of observables can be used to compute the SNR and Fisher matrix, which 
depend only on the available measurements; in [51 we give explicit expressions for the 
noise- orthonormal bases that simplify these computations; in [6] and [3 we apply our 
framework to the problem of determining the six-link, five-link, and four-link LISA 
sensitivity to sinusoidal signals and its precision (i.e., the expected statistical error 
due to instrument noise) in the determination of the polarization of monochromatic 
binaries, thus quantifying the heuristic statement that additional TDI observables help 
disentangle the GW polarization states. 



2. The LISA noise response 



We adopt the "classic" conceptualisation of the LISA measurement developed by 
Armstrong, Estabrook, and Tinto (see especially [3]), whereby: (i) each LISA 
spacecraft contains two optical benches, each with a proof mass and a laser; (ii) 
the inter-spacecraft one-way Doppler measurement "y" on each bench consists of the 
fractional frequency difference between the incoming laser (bounced off the local proof 
mass) and the local laser (unbounced); (iii) the intra-spacecraft measurement "z" on 
the same bench consists of the fractional frequency difference between the local laser 
(unbounced) and the laser from the other bench on the spacecraft, transmitted via 
optical fibre (and bounced on the other bench's proof mass). The noise in each of the 
measurements is modelled as due to the frequency fluctuations of the six lasers, to 
the random displacements of the six proof masses and of the six optical benches, and 
(for the y's) to shot noise and other optical-path noise in the low-SNR inter-spacecraft 
laser links. The optical-fibre noise can be removed by using the z's only in differences 
of the z's on the same bench; and the optical-bench motions along the lasers' lines of 
sight can be absorbed into the laser frequency noises. Without loss of generality, we 
assume also that the central frequencies of the lasers are all the same, and that all the 
Doppler beat notes due to the relative motions of the spacecraft have been removed 
via heterodyne measurements [6l [7] . 

Using the notation of [8] and [9], the y and z noise responses are given by 

--ft) = / " Li{t))-Cr(t) + y°l{t) - 2pmr{t) for unprimed /, 

^''^ ^' \ Cs{t-U (t)) -C; (<) + yZ {t) - 2pm; (t) for primed I, ^ ' 

and 

.noise f i it) ~Crit) + 2pm;it) for uuprimcd Z, 

^' I Cr{t)-C;{t) + 2p7nr{t) for primed ^' 
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where ysir denotes the inter-spacecraft measurement obtained for the laser incoming 
from spacecraft s into spacecraft r, and travelhng along link / (see figured] each link 
/ sits across from spacecraft I in the LISA triangle, and takes a prime if sir is an odd 
permutation of 123); furthermore, Zsir denotes the intra-spacecraft measurement on 
the same bench as ysir (this traditional notation is somewhat unfortunate, since the 
indices s and I refer to distant spacecraft that have no role in the local measurement). 
In these equations, the Cs and C* represent the six laser frequency noises, the prrir 
and pm* the six proof-mass velocity fluctuations along the lasers' lines of sight, the 
y°ir ^^-'^ inter-spacecraft optical-path noises, and L; is the one-way light time along 
link / (the asterisks over the C and pm do not denote complex conjugation, but rather 
different noise variables). 

TDI observables are linear combinations of several y and z, appropriately time 
delayed so that all instances of laser frequency noise {Cr and C*) cancel out. It is 
useful to introduce time-delay operators Vi such that for any measurement x{t), 

Vix{t) = x{t~Li{t)), VmVix{t)=x{t-Li{t)-Lra{t-Li{t))), ... (3) 

and so on, with the shorthand x-^i{t) = 'Dix{t), x-rni{t) = 'Dm'Dix{t), and so on. The 
time-delay operators commute only if the Li arc not functions of time, in which case 
the delays are usually denoted by subscripts set off by commas. For simplicity, in this 
paper we shall concern ourselves only with this case, which corresponds to "modified" 
TDI observables. ("Second-generation" TDI observables are necessary to remove laser 
noise completely because of the "breathing" of the LISA arms [TOllllj: however, when 
compared to the corresponding first-generation TDI observables, they engender only 
negligible corrections to the GW and secondary-noise responses, which are the building 
blocks of the SNR and Fisher matrix. Moreover, the second-generation observables 
can be approximated as finite-difference derivatives of their first-generation versions 
[TT| [9] , and this effect factors out in the computation of the SNR and Fisher matrix. 
For these reasons, the results of this paper remain valid for second-generation TDI.) 

The problem of finding combinations of the y's and z's that cancel the six laser 
noises, Cs and C*, can be reduced [ij to the problem of finding combinations of the 
new variable^ 

y'slr=yslr + ^{Zrl's,l- Zslr), (4) 

which cancel the three equivalent laser noises = (Cs + C*)/2. This reduction 
removes the z's from consideration!! and leads to the noise responses 

^,noisc^^^ _ I C'sA -Cr + vjr " '^P^r + pms,i - pm; for uuprimed I, 
1 C'^ i — C'j. + y°f^ — 2pm* -I- pm* i — pm^ for primed I. 

For terseness, in the rest of this paper we shall contract each of the index triples 312, 
123, 231, 21'3, 32'1, and 13'2 to its middle link index alone, as shown at the top of 
figure [11 so that yi = y^ir and zi = Zsir. 

X These substitutions are essentially equivalent, yet slightly different from those introduced in 
with the advantage of leading to more symmetric noise responses. 

§ It would seem that working with the y' it becomes hard to enforce the constraint that the z be 
used only in the same-spacecraft differences 2231 — ^Z2'lt ^312 — ^13/2, 2123 — ^21'3 (ori spacecraft 
1, 2 and 3 respectively), but this happens naturally in TDI combinations because of the way that 
the laser noises enter the y' noise responses. This is easy to see from the viewpoint of geometric 
TDI [9], where y' arrows that begin or end at the same spacecraft always appear in the sequences 

\'''V' V' (which involve the z difference at the spacecraft) or \ (which involves neither of the z's 
at the spacecraft). 
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y3i2=yi y2V3=yr 





unequal-arm relay monitor beacon 

Michelson 



Figure 1. Top panel: the LISA inter-spacecraft measurements. Eacii arrow 
represents tlie frequency-comparison measurement yi = y^i,. taken on the 
spacecraft at the tip of the arrow between the incoming laser (moving in the 
direction of the arrow and experiencing the light time L;) and the local laser. 
Unprimed (primed) link indices correspond to counter-clockwise (clockwise) laser 
propagation. Bottom panel: first-generation TDI observables as synthesized 
interferometers [9]. The four kinds of four-link observables are distinguished by 
the direction of the links converging on one spacecraft (in this case 1), which 
makes it a "centre," "relay," "monitor," or "beacon." 



The recently proposed [H] LISA architecture with a single proof mass on each 
spacecraft and with "strap-down" measurements (whereby the lasers are not bounced 
on the proof masses, but additional measurements are taken of the position of the 
latter with respect to the optical benches) would lead to analogous, but different noise 
responses. We expect the results of this paper to persist in the new architectures with 
slight modifications of SI but careful verification is certainly indicated once the LISA 
architecture is finalized. 



3. TDI observables from linear algebra and the LISA failure modes 



When the delay operators commute, they have an especially useful representation in 
the Fourier domain, 

Viiif) = Aiiif), P,„Ai(/) = A„AzJ(/), ... (6) 

where A/ = exp^'^'-''^' . This representation turns the search for laser-noise-cancelling 
combinations into a problem of linear algebra. The laser-noise content of the y'l 
measurements is described by the equation 

y' = ^eC', 



(7) 



where y' is the vector [y[,y2, y'3, y'y , y'2' , v's' 



a=[C[,C!2,C',f, and 



( 


-1 


Ai\ 


A2 





-1 


-1 


A3 








Ar 


-1 


-1 





A2- 


^3' 


-1 


) 



(8) 
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now, a generic observable (not necessarily laser-noise-cancelling) is a linear 
combination a-^y' (where denotes vector transposition), where the coefhcients 
a will normally be polynomials of the A; with integer coefficients. It follows from ([7]) 
that the a that correspond to laser-noise-cancelling observables must span the null 
space of Dq, 

a^DeC' = Q ^ D^a = 0; (9) 

thus, the task of finding TDI observables reduces to finding a basis for null(r'l'), and all 
TDI observables can then be represented as linear combinations of the basis elements 
(although such representations may involve Aj's, and therefore delays, as we shall see 
in section 15. 2p . The Fourier-domain representation is especially suited to our goals of 
assessing the reduction in SNR associated with particular failure configurations, and 
of evaluating the parameter-estimation precision for observed gravitational signals, 
since these operations involve calculations performed in the Fourier domain. 

The loss of one or more inter-spacecraft measurements results in the loss of the 
corresponding y'^. However, the linear-algebra problem of (O can still be formulated 
after removing the appropriate rows in y' and Dq. In general, if n is the number of 
measurements and m the number of laser noises to be cancelled, and if the subsetted 
Dq has full rank (which is generally the case), it follows that the dimension of null(-D"^) 
is r = n — m. Thus, for "full LISA," with six inter-spacecraft measurements and three 
noises, we can construct three independent TDI observable^; these are reduced to two 
after one measurement is lost, and to one after two are lost. In section [5] we describe 
which of the standard TDI observables can be reconstructed in each case. 

The effects of losing one or more intra-spacecraft measurements are subtler. 
Consider for instance losing Z231 = 2:3 (on spacecraft 1): from ([4]), we see that we 
immediately lose 2/3 and j/g, ; we also lose the ability of using y'2 and y'2,, which 
involve the other measurement on spacecraft 1, Z32'i = 22' (remember that Z3 and 
Z2' can appear only in the combination Z3 — z^i if fibre noise must be cancelled; thus, 
losing one or two z measurements on the same spacecraft is equivalent). However, it 
turns out that the combinations y'^, -\- y'^ 3/ and j/j + 1/2' 2 graphic notation of 

geometric TDI and p') contain neither Z3 nor Z2' (and also no C[i}. Thus we 

have four available measurements and two noises to be cancelled, so two independent 
TDI observables can be constructed. If some or all z measurements are lost on two 
spacecraft (say, 2 and 3) then only two combined measurements are available (in this 
case, j/3 + j/g, 3 and y'21 +2/22') with one laser noise to cancel, so a single independent 
TDI observable can be constructed (it is the unequal-arm Michelson centred on the 
non-faulty spacecraft). If both inter- and intra-spacecraft measurements are lost, then 
at most one TDI observable can be constructed (and it is always an unequal-arm 
Michelson), and in some cases none. Again, see section [5l It is also possible to come 
to the same conclusions by working with a larger linear system that involves the 
available yi (not ?/j), zi and the six laser noises Cr and C*. 

4. Invariance of SNR and Fisher matrix by change of TDI basis 

Of course, the basis of null(D-'") (with D the appropriate row-reduced matrix for the 
available measurements) is not unique; however, we now show that the total SNR 

II Although a combination of the three is insensitive to GWs in the long-wavelength limit because of 
symmetry; sec figure [2] 



Sensitivity and parameter- estimation precision for alternate LISA configurations 6 



and Fisher matrix are invariant with respect to the choice of a TDI basis, so they 
depend only on which yi and zi measurements are available, and can computed over 
any convenient basis. From the SVD decomposition '1 3] 

= ^ ■ ^ ■ (with S.y = for i 7^ j), (10) 

mxn mxm mxn nxn 

it foUows that a basis for nuU(Z3-'") is given by the columns of V that correspond to 
the "extra" r null columns of S. The SVD decomposition is degenerate with respect 
to nonsingular linear combinations of these columns (i.e., with respect to choosing 
a different TDI basis). Calling W the n x r matrix of these columns, the optimal 
matched-filtering SNR for a source h(/) whose components are given by the individual 
TDI observables Xa — yi{t)Wia is given by 

r+oo 

SNR2 = 4Re / hi{f)Sj{f)hpif) df (11) 

JQ 

r+oo 

= 4Re/ h^if)[W^Cif)WT'Hf)df, 
Jo 

since 

lso.p{f)6{f - /') = mf)w,o.][y',{f')w,pr) = wM{f)y']{n)w;0 (12) 

= \w^C{f)W*5{f-f'), 

where the dj components are the secondary-noise cross-spectra of the individual y'^ . 
As discussed in [Appendix A[ Cij matrix is in general complex and Hermitian. 

If we change the TDI basis by means of a nonsingular linear transformation 
h —^ Bh (so that we are computing Sa'p' over the TDI observables X'^ = Ba'aXa — 

Ba'ay'i{t)Wia) WC get 

h''[W^CW*]-^h h''B''[BW^CW*B^]-'^Bh = B\B^)-^[W'^CW*]-^B-^Bh 

^h''[W'^CW*]-^h, (13) 

so the SNR is invariant. Indeed, any noise inner product (h, g), and therefore the 
Fisher matrix elements (h_^,h^i,), are also invariant under similar transformations. 

As a corollary we get the intuitive fact that computing the SNR with an extra 
TDI observable that is a linear combination of the r others (which amounts to a 
(r + 1) X r matrix B) does not change the result. To see this, we inject an additional 
(r-f 1) X (r + 1) linear transformation B' that makes one of the (r-|-l) TDI observables 
identically null; the entire product (jl3p can then be rewritten by dropping one of the 
rows of B' B, so we are back to the case of a nonsingular r x r transformation. 

The equivalence of different TDI bases for the purpose of computing the Fisher 
matrix can also be proved with a different approach, closer to the reasoning of Romano 
and Woan [14] : we consider the estimation problem of determining the GW parameters 
and the laser frequency noises in the presence of secondary instrument noise, assumed 
Gaussian and stationary. The likelihood of the LISA data is then written in terms 
of the y- and their secondary-noise cross-spectra, and (without loss of generality) 
rewritten in a TDI/non-TDI basis where the first m components span mil\{D^) and 
the other r span null(I?^)^. We compute the Fisher matrix by taking derivatives of 
the log Ukelihood with respect to the GW parameters and to the individual frequency 
components of laser noise (again, our premise is that we are estimating these in 
addition to the GW parameters, even if in practice it will turn out that only some 
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combinations of the laser noises can truly be estimated from the LISA measurements). 
Using the Frobenius-Schur formula [T5], we can then show that the GW-parameter 
sector of the inverse Fisher matrix (i.e., the expected covariance matrix for the GW- 
parameter estimators) is equal to the inverse of the TDI sector of the full Fisher 
matrix, with all the laser-noise components dropping out. Since the projection on the 
null(D"^) and null(D"^)-'- subspaces can be cast as a geometrical operation, we see that 
the Fisher matrix does not depend on the choice of TDI observables, but only on the 
geometry of null(Z?"^). 

5. Noise-orthogonal bases of TDI observables 

An especially useful linear transformation of a TDI basis is the one that leads to an 
orthonormal basis of Sap, since the SNR then has the simplified form 

SNR2 = E 4 Re ^iZl^d/; (14) 

^, Jo Ja'a'[J) 

a similar expression follows for the Fisher-matrix elements. We now work out such 
bases and the corresponding Sa'f!' under various LISA failure modes. 

5.1. Four-link configurations 

We begin with the case of two lost yi measurements, where null(Z)"^) has dimension 
one, so that essentially only one TDI observable can be constructed, which is trivially 
a noise-orthonormal basis. Depending on which measurements are lost, this observable 
is one of the standard four-link observables (unequal-arm-Michelson, relay, beacon, 
and monitor) described in [TT], and shown pictorially in the bottom part of figure 
[1] In this section we spell them out and identify them with a new compact naming 
scheme. 

If the missing yi are along the same arm (say, yi and y'l), we can build the 
unequal- arm Michelson observable centred on spacecraft 1, 

Xl = 2/2' + 2/2,2' + 2/3,22' + 2/3',322' " iVi + V's' ,3 + 2/2',3'3 + 2/2,2'3'3); (15) 

the analogous X2 and X3 are obtained by the cyclic permutations 1-^2^3 and 
1' ^ 2' ^ 3' (which will go without saying for all observables to follow). Note that 
these Xi, X2 and X3 are usually called X, Y and Z in the TDI literature, but here 
we prefer a notation that emphasizes which yi enter them. 

If the missing yi are consecutive and "codirected" (say, 2/2' and yy), we can build 
the "forward" relay observable Ui that goes through spacecraft 1 by way of links 3 
and 2, 

Ul = y'v + V'l.v + 2/24I' + 2/3,21'l - (2/2 + 2/3,2 + 2/i,32 + 2/'i',i32); (16) 

by contrast, the "backward" relay observable Vi goes through spacecraft 1 by way of 
2' and 3', 

^1 = 2/'l + 2/i',l + 2/3', I'l + 2/2',l'13' - (2/3' + 2/2',3' + 2/i',2'3' + 2/'l,l'2'3')- (1^) 

If the missing yi are consecutive, but have opposite "directions" , we have two 
cases. If the available measurements are (say) y\ and yy , and also 2/2' and 2/3 (which 
"point" towards spacecraft 1), we can build the monitor observable 

El=y'3+ 2/1,3 + 2/l',31 - 2/2' - 2/l',2' - 2/i,2'l' + 2/2', I'l - 2/3,i'i; (18) 
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Table 1. Standard TDI observables available in various LISA failure modes. For 
a single lost yi , all TDI observables on the corresponding line (or column) of the 
left table are available (and any two of them are independent); for two lost yi, 
only the observable at the corresponding intersection is available. For one or two 
lost zi on the same spacecraft, the two TDI observables on the corresponding line 
(or column) of the right table is available; for two to four lost z; on two spacecraft, 
only the observable at the intersection is available. For mixed failure modes, at 
most one Xi observable is available, as can be seen by intersecting the two tables. 



yi y2 ys yv y2' ys' \ Zi on i z; on 2 z; on 3 

Vs V2 Xi Es Pa Zi on 1 X3 X2 

V3 Vi Pi X2 El zi on 2 X3 Xi 

V2 Vi E2 Pi X3 zi on 3 X2 Xi 

Xi P3 E2 U3 U2 

E3 X2 Pi U3 Ui 

P2 El X3 U2 Ui 



yi 
2/2 
2/3 

2/1' 
2/2' 
2/3' 



whereas if the available measurements other than yi and yy are 2/2 and ?/3' (which 
"point" away from spacecraft 1), we can build the beacon observable 

-Pi = 2/1,2 + yi',21 + y3',211' - y'l',3' - 2/1,3'!' - 2/2,3'l'l + 2/2,3' " 2/3',2- (19) 

Again, the analogous observables i?2, E-^, P2, and P3 can be obtained by cyclic 
permutations. The left part of table [1] shows which observable is available in each 
case: look at the intersection of the row and column corresponding to the missing yi . 

We can now compute the noise PSDs for these observables by writing them in the 
Fourier domain in terms of the yi and of the complex A; , multiplying them by their 
complex conjugates, and substituting the noise responses ([5]). We work in the limit 
of equal armlengths (so all A; — exp 2n\fL) and assume that all noises are Gaussian 
and uncorrelated (i.e., have null cross-spectra), and have spectral densities 5'°p(/) for 
the optical-path noises and 5p™(/) for the proof-mass noisesl^ The resulting PSDs 
are 

Sx, if) = 16 5°P sin^ x + 16 SP'^iS + cos 2x] sin^ x, (20a) 
Su, (/) = 5y, (/) = 8 S°P [4 + 4 cos x + cos 2x] sin^ f 

+ 16S'P'"[5-h5cosa;-t-2cos2x]sin2 |, (205) 
SEiif) =Sp,{f) = 8S'°P[3 + 2cosa;]sin2f + 16S'P'"[3-Hcosa;]sin2f, (20c) 
where x — 2TrfL. 

5.2. Five-link configurations 

If one yi is lost, null(I?"^) has dimension two, corresponding to two independent TDI 
observables. Depending on which yi is missing, a different set of the five standard 
observables described in the last section can be constructed, as given by the rows of 
the left part of table [1] Within each set, any two observables can form a basis for 
null(-D"^). This means also that any two observables can be used to re-express any 
other. To find such relations, we solve linear systems such as 

Xi=c{Xi,V2)V2+c{Xi,V3)V3, (21) 

^ For reference, the standard LISA model used in the Mock LISA Data Challenges 1161 has 
(S°P)i/2 = 20 X 10-12 mHz-i/2_ (5Pm)i/2 ^ 3x io-i5[i + (io-'iHz//)2]i/2ms-2Hz-i/2 and 
L = 16.6782 s. 
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and therefore 

/ \ 



1- A3, A3 
-1 + A2A2' 

\-A3 + A3A2A2' / 

which yields 



/ 

-1 + A2'A2 - 

1- A2'A3,Ar 
-A3/A1/ + A2 


V-Ai- + Ai-A2'A2 



\ 

-A2' + A2'A3'A3 



-1 + A3-A3 
1- Ai,A2'A3, 
-Ai,A2' + A3 



(22) 



X, = 



A2' - A2'A3'A5 



1 - A3'Ai'A2 



A2A2 



1 



1- A3'ArA2 



-v. 



and finally 



Analogously, 



Xi - Xi, 



3'1'2' 



V2. 



2,2'3'3 



V3. 



V3. 



^1,132 — Xi = C/2,3'3 — U2 + C/3,3 ^ C^3,322', 
-'^'1,32 — Xi^i' = P2 — -P2,3'3 + £^3,322' — £'3,3, 
-'^1,3'2' — -'^l,! — -E'2,2'3'3 — E2.2' + P3 ^ £3, 22' 



(23) 

(24) 

(25 a) 
(256) 
(25 c) 



and also 



^2,1' — V2,32 — £2,3'!' — ^2,2 + £3,1' ^ £'3,1'2'2, (26a) 
^34' ~ ^3,32 = P2 ^ -f2,3'3 + -E'3,1'2' ^ £3,3, (26&) 
C/2,2'3' - U2,l = £2,31 - E2,2' +P3- ^3.22', (26c) 

C/3,2'3' — C^3,i = E2.1 — £2,133' + £3,21 — £3,3'- (26rf) 

(Since these equations cease to hold strictly if the delays become noncommutative, 
here the ordering of the delay indices has little actual meaning.) 

We now seek a basis of two noise-orthogonal TDI observables by diagonalizing 
the cross-spectrum matrix Saf3 of two of the five standard observables, under the 
assumptions on the noises given above. The diagonalisation process has different (but 
equally valid) results depending on which two we use, but not all choices are equally 
convenient. For instance, if t/i is lost and we work with P2 and £3, we obtain a 
correlation matrix that has two equal terms on the diagonal [see (j20cp ] and complex 
cross terms 



Sp^E, =S*E,p, =4sinf sina;e'^/2[5.op(3_^gix-)_^25P'"(3-e-'^)]. (27) 

Unfortunately, the resulting eigenvectors are oc £3 — &'^'^^^P2 and £3 -I- &'^^^'>P2, with 
0(/) — arg^Sp^Es)', that is, they have complex and frequency- dependent coefficients, 
whereas we are used to orthogonal observables that are real variables and have integer 
coefficients. The same problem occurs if we diagonalize the cross-spectrum matrix for 
any two out of Xi, V2, V3, P2, £3. After some experimenting, we find that the linear 
combinations P2 — Xi/2 and £3 -I- A'i/2 do have a real cross-spectrum; diagonalizing 
them yields the eigenvectors 

Ad) - (F2 

with PSDs 



E3~Xi)/V2, £(1) = (P2+£3)/V2, 



(28) 



Sa(i) =8 5°P[2-f cosxl^sin^^ 



Sew 



8S'°P[2- 



cos a; sm 



2 

2 X 

2 



8 5P''^[14 
32S'P'"[2- 



15 cos X -f 6 cos 2a; -I- cos 3a;] sin^ |, (29a) 
cosa;]2sin2|. (296) 
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Working through a similar exercise when yii is lost yields 

^(l') = (P3-i?2-^l)/V2, = (P3+i?2)/%/2, (30) 

with the same PSDs. Once again, cyclic index permutations will yield orthonornial 
basis for the other missing yi. 

5.3. Full six-link configuration 

Here null(i?|") has dimension three, and all standard observables can be constructed. 
An obvious choice is to diagonalizc the three unequal-arm Michelson observables, 
which have cross-spectra 

Sx^x, = -8 5°P cos X sin^ x - 32 S^"^ cos x sin^ x [i ^ j) (31) 

and yield the well-known eigenvectors 

A(^) = (X3 - Xi)/x/2, E^^^ = (Xi - 2X2 + X3)/\/6, 

tW = (Xi+X2 + X3)/V3, (32) 

with PSDs 

Saw =Se(x) = 8S'°p[2 -f cosx] sin^x 16S'P"'[3 2cosx cos2a;] sin^x, (33a) 
Stix) = 16 S'°P[1 - cosx] sin^ x + 128 S^p™ sin^ x sin'' f . (336) 

Less well known is the fact that Ei, E2, and E^ , and Pi, P2, and P3 are also suitable 
bases. Since 

Sp^P^ = Se,e, = 4 ^"P sin^ X - 8 SP'''[-1 + 2 cosx] sin^ x, {i ^ j), (34) 
the resulting eigenvectors have structure similar to ((32|) . 

A(^) = {E3 - Ei)/V2, P'^) = {El - 2E2 + E3)/V6, 

T^^^ ^(Ei+E2+E3)/V3, (35) 



and 



A(^) = (P3 - Pi)/V2, = (Pi - 2P2 + P3)/V6, 

t(^) = (Pi + P2 + P3)/V3, (36) 

with PSDs 

Sa(e.p, = Se(b,p) = 85°p[2 + cosx] sin^ | + 16S'P'"[3 -|- 2cosx -t- cos2x] sin^ f , (37a) 

St(e,p) =8S'°P[5 + 4cosx]sin2|-f 32S'P'"[5 + 4cosx]sin''f. (37&) 

By contrast, working with Ui, U2, and U3 (or Vi, V2, and V3) leads to complex 
cross-spectra (although rather symmetric, with SuiU2 = SU2U3 — SU3U1 = ^U2Ui ~ 
^U3U2 ~ ^UiUa)^ ^^'^ ^^'^ '^^^ three of the resulting eigenvectors have complex 
coefficients, and do not correspond to real observables (the third eigenvectors are the 
again completely symmetric {Ui + U2 + C/s)/^^ and {Vi -I- V2 -I- Vi3)/\/3). 

The classic noise-orthogonal observables are perhaps the A, E and T of \17\ I18j . 
written in terms of the first-generation-TDI Sagnac observables 

a = y'2' + y'v.2' + y-A\vi' - {y's + yi.3 + (38) 

(with /3 and 7 obtained by cyclical permutations) which are a basis for null(£)"^) only 
if A; = A;/ , and have PSDs 

Saa = Sfifj = S-^^ = 6 S'°P + 8 S'P"' [5 4 cos X -I- 2 cos 2x] sin^ f , (39a) 
S^p = Sf3j = Sja^2 S°P [2 cos X cos 2x] 4 ^p'" [cos x - 1] , (39 &) 
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yielding the eigenvectors 

A^{j-a)/V2, ^ = (a - 2/3 + 7)/V6, 

T = {a + l3 + l)/V3 (40) 

with PSDs 

Sa = Se = 8S'°P[2 + cosa;]sin2 f + 16 S'p™[3 + 2 cosx + cos 2a;] sin^ |, (41a) 
St = 2S'°p[1 + 2cosa;]2 + g^P^'sin^ ^. (416) 
If A; 7^ All we need the more compHcated 

ai — Z/2' + 2/l',2' + 2/3', 1'2' + y3,3'l'2' + yi,33'l'2' + J/2,133'1'2' 

- + yi,3 + ^2,13 + y2',213 + 2/l',2'213 + y3M'2'213) (42) 

(and similarly for a2 and a^), which again have (I32p -like eigenvectors A, E and T, 
and PSDs 

S*^ = 4sin2^ X 5*^, = 4sin2f^ X Sb, (43) 
since in the limit A; A;/ ai ~ a — ai23, and so on. 

5.4- Missing-zi configurations 

For lost zi on two spacecraft we go back to the unequal-arm Michelson four-link 
scenario, as we do for certain combinations of lost yi and zi] the only scenario that 
remains to be covered is one of lost z/ on a single spacecraft (say, 1), where two 
unequal-arm Michelson observables (in this case, X2 and X3) can be constructed. 
From (|20ap and (pij) we get the eigenvectors 

= {X2 - X3)/ V2, £;(^i) = {X2 + X3)/V2, (44) 

and the PSDs 

SA(.i)^Ei.^) = 85°P[2±cosx]sin2x + 85P'"[6±4cosa; + 2cos2a;]sin2a;; (45) 
in the ± signs, the + refer to A'^-'^', the — to E'^'^^\ 

6. The LISA sensitivity v^rith four, five and six inter-spacecraft 
measurements 

To illustrate the use of noise-orthogonal TDI bases, we now compute the sky-averaged 
LISA sensitivity to monochromatic, sinusoidal signals using four, five and six yi 
measurements. We adopt the expressions of Vallisneri [8j for the LISA orbits and 
its response to GWs, and we consider the GWs emitted by circular, nonspinning 
binaries, 

h+{t) = h{l + cos^ t) cos(27r/< + if a) 

/ix(0 = -2/1 cost sin(27r/< + (po), (46) 

where l is the inclination of the binary's orbital plane with respect to the line of sight 
to the Solar-system barycentre (SSB), / is the frequency of the source as measured at 
the SSB, (fo is the initial phase and h is the overall GW strength. 

The sensitivity to sinusoidal signals, for sources at a fixed position in the sky, 
has been defined traditionally as the GW strength h "required to achieve a SNR of 
5 in a one- year integration time, as a function of Fourier frequency" [in] ■ We follow 
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|19| in computing not the sky-average of this position-wise sensitivity, but instead the 
GW strength h required to achieve an rms-averaged SNR of 5 for sources isotropically 
distributed over the sky. However, unhke [12], where polarization states are chosen 
uniformly over the abstract Poincare sphere [20] . we consider sources isotropically 
distributed over t and over the GW polarization ip (essentially a rotation of /i+ and 
hx with respect to the principal axes conventionally defined at every sky position). In 
fact, we find empirically that the two distributions yield almost identical sky-averaged 
sensitivity curves. 

The rms sky average has the advantage that the LISA orbital motion can 
be neglected, since the source distribution remains isotropic (and the sky average 
invariant) as the LISA orientation changes along the year; thus, all SNRs can be 
computed for a stationary LISA sitting at the SSB. This approximation neglects the 
Doppler modulation of individual sources, which changes results only marginally by 
smoothing out the signals across nearby frequencies. Following [T9 j , we can then 
operate in the frequency domain, as follows. For 10, 000 uniformly sampled sky 
positions, inclinations and polarization states, we turn h^{t) and hx{t) into complex 
phasorO that are inserted in the yi GW responses (l)-(2) of [8^, which in turn are 
inserted in the TDI expressions for the observables; delays are always replaced with 
A; factors, and evaluated in the limit of equal armlengths. The squared modulus of 
the resulting phasor (which is a function of frequency through the A;), multiplied 
by T = 1 year and divided by the noise PSD, yields the SNR^, which can then be 
summed over noise-orthogonal observables and averaged over sky position, inclination 
and polarization. 

Our results are shown in figure [2l As plotted in the bottom panel, the four- 
link configurations yield slightly different sensitivities depending on which type of 
observable can be constructed (Xj; Ui and Vi, which have the same sensitivity; or Ei 
and Pi, which also have the same sensitivity); however, they all converge in the low- 
frequency limit, and they essentially agree at high frequencies. The gain in upgrading 
to a five-link configuration is \/2 in the low- frequency limit, and as much as ^ -y/S 
at peak frequencies above 10 mHz. The further gain in upgrading from a four-link 
to a six-link configuration is again \/2 in the low-frequency limit (because the third 
noise-orthogonal observable T is insensitive to GWs in that limit, as seen in the top 
panel), and as much as ~ 2 at peak frequencies above 10 mHz. 

7. Polarization-angle estimation error with four, five and six 
inter-spacecraft measurements 

As a second example of the use of noise-orthogonal TDI bases to characterize the 
LISA performance under different failure scenarios, we now compute the expected 
estimation error for the GW polarization ip of our fiducial monochromatic binaries, as 
predicted by the appropriate diagonal element of the inverse Fisher matrix (see, e.g., 
[21]). One of these binaries is completely described by 7 parameters (ecliptic latitude 
and longitude, polarization, amplitude, inclination, frequency, and initial phase) so the 
full Fisher matrix that we invert is 7 x 7. To compute its elements, we work in the time 
domain with the full LISA GW response of [8| (including the amplitude and Doppler 
modulations due to the LISA orbital motion), and we compute signal derivatives with 
respect to source parameters by means of finite differences for very small parameter 

+ With amplitudes reduced by 1/2 to account for a constraint of reality on the GWs. 
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Figure 2. Bottom panel: LISA sensitivity to monochromatic sinusoidal signals 
in the six-link, five-link, and in all four-link configurations (which are different for 
Xi and for the two Ei/Pi and Ui/Vi pairs). Top panel; sensitivity of the individual 
noise-orthogonal observables in the six- link E(^) and Tt^') and five-link 

configurations, as compared to the standard unequal-arm Michelson 
observable Xi. 

displacements. When more than one noise-orthogonal observable can be constructed, 
the corresponding Fisher matrices are summed before inverting to yield errors. 

We consider 10, 000 binary systems spread across the LISA band, with uniformly 
sampled sky positions, inclinations, and polarizations; furthermore, we consider 
observations lasting one month, three months and one year, with four- link (^^^i), five- 
link + E^^^) and six-hnk (A^-^' + E^^'^ + T^^)) LISA configurations. Figure [3] 
shows the median polarization errors as a function of source frequency. The signal 
amplitudes were chosen individually for each system so that it would be detected with 
an SNR of ten by the full six-link LISA configuration, notwithstanding the observation 
timespan (thus, the one-year signals are intrinsically weaker than the three-month 
signals, and these in turn are weaker than the one- month signals). Several features 
are worth discussing: 

• Errors improve in the longer observations, because the LISA orbits (which have 
a period of one year) can then modulate the LISA response more strongly, 
helping to disentangle the sky-position, polarization and inclination source 
parameters. This improvement is most noticeable between one-month and three- 
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Figure 3. Median polarization uncertainty as a function of binary frequency, 
for 10, 000 sources uniformly distributed over sky positions, inclinations, and 
polarizations, as determined using four-, five- and six-link LISA configurations. 
The three panels show expected uncertainties for observations of 1 month, 3 
months and 1 year. The signals are scaled so that their six-link SNR is 10; 
thus, the four- and five-link uncertainties are computed for signals that would 
have four- and five-link SNR less than 10. 

month observations, and especially so for the four-link configuration, for which 
the one-month errors are so large (~ tt/2) that they are off the chart in this figure, 
and that polarization is essentially undetermined at SNR = 10. 

• The five- and six-link errors coincide below 1 mHz, just as the five- and six- 
link SNRs do. Furthermore, above 1 mHz the ratio of the six-link to five-link 
errors agrees (within numerical noise) with the ratio of the six-link to five-link 
SNRs. Recalling that Fisher-matrix expected errors scale as SNR~^, we see that 
the primary effect of switching from five- to six-link configurations is to improve 
estimation error by increasing the SNR, but not by providing complementary 
"views" of the same signal with different geometries. (Thus, if the five-link curve 
is renormalized by setting the five-link SNR of all binaries to 10, it collapses onto 
the six-link curve.) Interestingly, the SNR-renormalized five- and six- link errors 
are almost constant between 0.1 and 10 mHz. This fact was already noticed by 
Crowder and Cornish [22] ; unfortunately, explaining how it comes about requires 
delving into the analytic structure of the LISA response, and is outside the scope 
of this paper. 

• The four-link errors behave a bit differently: their ratio to the five- and six-link 
errors is also explained by the SNR~^ scaling, but only for observations of three 
months or longer. Below three months, the four-link errors get worse and worse 
(relative to the five- and six-link errors) as the observation time gets shorter. 
We observe this behaviour in figure [H which plots the SNR-renormalized four- 
and five-link errors for binaries with GW frequencies between 1 and 1.1 mHz, as 
functions of observation time. The curves merge at three months, but the four- 
link errors degrade much more quickly than the five-link errors as observations get 
shorter. The apparent reason is that when we reduce the observation time, the 
polarization and sky-position parameters become more degenerate with four links 




Sensitivity and parameter- estimation precision for alternate LISA configurations 15 



1.4 




I ' ' ^^—^^^—^—i 

1.0 1.5 2.0 2.5 3.0 

observation time [months] 



Figure 4. Median polarization uncertainty as a function of observation length, 
for 10, 000 sources uniformly distributed over sky positions, inclinations, and 
polarizations, as determined using four- and five-link LISA configurations. The 
signals are scaled separately before computing the four- and five-link errors so 
that their four- and five-link SNRs are both 10. 

than with five; correspondingly, the parameter errors become more correlated, 
and therefore individually larger. Thus the additional "view" of the signal that 
is available with five links does help the polarization error, but only for rather 
short observations; above three months, its effect gets washed out by the greater 
amount of information provided by the longer orbital baseline. 

8. Conclusion 

We have discussed a simple framework to assess the LISA science performance (as 
indexed by its sensitivity to sinusoidal signals and by the inverse-Fisher-matrix 
parameter errors) in its full six-link configuration, and in failure modes where up 
to two inter-spacecraft laser measurements and up to four intra-spacecraft laser 
measurements become unavailable. We have described a unified, straightforward 
procedure to find the TDI observables that can be constructed in each failure 
mode, and to derive noise-orthonormal bases for them. In particular, we have given 
explicit expressions for the ^f*-' and i?'*-* orthonormal observables possible with five 
links, which had not previously appeared in the literature. Furthermore, we have 
demonstrated that both the SNR and the Fisher matrix are invariant under changes 
of TDI bases, so they depend only on which measurements are available, and can be 
computed for any convenient choice of observables. As an example of the application 
of our framework, we have computed the LISA sensitivity to sinusoidal signals in 
six-link, five-link, and four-link configurations (figure ^ and the expected precision 
for the estimation of the polarization angle of monochromatic binaries in the three 
configurations (figure [3]) . The framework assumes a "classic" conceptualisation of 
the LISA noises and measurement, but we expect that our results could be adapted 
easily to the emerging "strap-down" architecture without significant changes. The 
framework can also be extended trivially to derivative LISA missions that include 
more than three spacecraft (such as the proposed "bow-tie" LISA [53]) and possibly 
more than six inter-spacecraft measurements. 
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Appendix A. Multi-observable formulation of the noise inner product 

The generalisation of the standard noise inner product to multivariate processes 
(in this case, multiple TDI observables) follows from writing the probability of the 
stationary Gaussian vector process n„ as 

pih) (X exp h*MY^^S-'^{f)hpif)df^ EE exp{-(h,h)/2}. (A.l) 

which is correct since the correlations are entirely described by the cross spectrum 
{'^a{f)n*ij{f')) — '•^^Sai3{f)S{f — /'). Using the fact that Sap must be Hermitian, and 
furthermore that Sapi^f) = S^/^if) (since n{—f) = n*{f) for real noise processes), 
and replacing '^'^^Saj3{f) with the more familiar Sap = '^'^^^Sap (which becomes the 
one-sided spectrum for a = f3), we can write 

r+co 

(g,h)=4Re/ 9:{f)S-'pif)hi3{f)df. (A.2) 
Jo 

Note that the imaginary part of Sa/s is crucial to computing this inner product 
correctly. Saf3 can only be real if the correlation function Cap (t) = "^a {t + 

T)np{t)dt is even, but this is never the case for TDI variables, which include time- 
delayed linear combinations of different noises. Consider for instance two noise 
processes ni{t) and n2{t) = ni{t — At): 

ni{t + T)ni{t-At)dt = / ni(t + r + Ai)ni(t) dt = Cii(r + Ai).(A.3) 

-oc J — CO 

Since for a stationary process rii the autocorrelation Cii(t) is even, Ci2(t) cannot 
be even. It is true however that the imaginary part of Sap may be disregarded 
in computing the cross power (i.e., integrating Sap over all frequencies) for real 
processes, since then all imaginary contributions cancel between positive and negative 
frequencies. 

Note that in this paper we adopt the Numerical Recipes [24] definition of the 
Fourier transforms, namely h{f) — J h{t) exp^^'-''* dt, h{t) = J h{f) exp^^^'-'^* dt. 
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